Reciprocal expression of MADS-box genes and DNA methylation reconfiguration initiate bisexual cones in spruce.

The naturally occurring bisexual cone of gymnosperms has long been considered a possible intermediate stage in the origin of flowers, but the mechanisms governing bisexual cone formation remain largely elusive. Here, we employed transcriptomic and DNA methylomic analyses, together with hormone measurement, to investigate the molecular mechanisms underlying bisexual cone development in the conifer Picea crassifolia. Our study reveals a "bisexual" expression profile in bisexual cones, especially in expression patterns of B-class, C-class and LEAFY genes, supporting the out of male model. GGM7 could be essential for initiating bisexual cones. DNA methylation reconfiguration in bisexual cones affects the expression of key genes in cone development, including PcDAL12, PcDAL10, PcNEEDLY, and PcHDG5. Auxin likely plays an important role in the development of female structures of bisexual cones. This study unveils the potential mechanisms responsible for bisexual cone formation in conifers and may shed light on the evolution of bisexuality.

Diploid species phylogeny and evolutionary reticulation indicate early radiation of Ephedra in the Tethys coast.

Reconstructing a robust species phylogeny and disentangling the evolutionary and biogeographic history of the gymnosperm genus Ephedra, which has a large genome and rich polyploids, remain a big challenge. Here we reconstructed a transcriptome-based phylogeny of 19 diploid Ephedra species, and explored evolutionary reticulations in this genus represented by 50 diploid and polyploid species, using four low-copy nuclear and nine plastid genes. The diploid species phylogeny indicates that the Mediterranean species diverged first, and the remaining species split into three clades, including the American species (Clade A), E. rhytidosperma, and all other Asian species (Clade B). The single-gene trees placed E. rhytidosperma sister to Clade A, Clade B, or Clades A + B in similar proportions, suggesting that radiation and gene flow likely occurred in the early evolution of Ephedra. In addition, reticulate evolution occurred not only among the deep nodes, but also in the recently evolved South American species, which further caused difficulty in phylogenetic reconstruction. Moreover, we found that allopolyploid speciation was pervasive in Ephedra. Our study also suggests that Ephedra very likely originated in the Tethys coast during the late Cretaceous, and the South American Ephedra species have a single origin by dispersal from Mexico or North America.

Tracing the geographic origin of endangered plant species using transcriptome-derived SNPs: An example of Cathaya argyrophylla.

Genetic markers have emerged as one of the most promising tools for species identification and geographic traceability in biodiversity conservation and international trade of biological products. However, traditional molecular markers rarely have sufficient resolution at lower taxonomic levels, especially for discriminating closely related forest tree species and their populations. In this study, we developed a panel of RNA-Seq based single nucleotide polymorphism (SNP) markers for tracing the geographic origin of an endangered conifer, Cathaya argyrophylla, which is a paleoendemic restricted to four mountain regions in subtropical China. A total of 69 individuals from five populations (DLS, SHS, HP, BMS, and DYS) covering the entire range were used for transcriptome sequencing. Based on these transcriptomic data, we evaluated genetic variation and population structure of C. argyrophylla, and found extremely low nucleotide diversity but strong population differentiation. We also screened 113 population-specific SNP loci, including 96 for BMS, eight for DYS, six for SHS, two for HP, and one for one of the three subpopulations from DLS. According to these geographically diagnostic SNPs, we designed four population-specific molecular barcodes for PCR amplification. To test the utility and efficiency of the four markers in geographic discrimination, double-blind experiment was performed using 157 individuals labelled without any locality information. We found that almost all tested individuals could be successfully assigned to their geographic localities. Our study not only sheds some new light on the genetic profile of C. argyrophylla, but also provides a practical and cost-efficient solution for geographic traceability using transcriptome-derived SNPs.

A Comprehensive Evolutionary Study of Chloroplast RNA Editing in Gymnosperms: A Novel Type of G-to-A RNA Editing Is Common in Gymnosperms.

Although more than 9100 plant plastomes have been sequenced, RNA editing sites of the whole plastome have been experimentally verified in only approximately 21 species, which seriously hampers the comprehensive evolutionary study of chloroplast RNA editing. We investigated the evolutionary pattern of chloroplast RNA editing sites in 19 species from all 13 families of gymnosperms based on a combination of genomic and transcriptomic data. We found that the chloroplast C-to-U RNA editing sites of gymnosperms shared many common characteristics with those of other land plants, but also exhibited many unique characteristics. In contrast to that noted in angiosperms, the density of RNA editing sites in ndh genes was not the highest in the sampled gymnosperms, and both loss and gain events at editing sites occurred frequently during the evolution of gymnosperms. In addition, GC content and plastomic size were positively correlated with the number of chloroplast RNA editing sites in gymnosperms, suggesting that the increase in GC content could provide more materials for RNA editing and facilitate the evolution of RNA editing in land plants or vice versa. Interestingly, novel G-to-A RNA editing events were commonly found in all sampled gymnosperm species, and G-to-A RNA editing exhibits many different characteristics from C-to-U RNA editing in gymnosperms. This study revealed a comprehensive evolutionary scenario for chloroplast RNA editing sites in gymnosperms, and reported that a novel type of G-to-A RNA editing is prevalent in gymnosperms.

Phylogeny and evolution of Cupressaceae: Updates on intergeneric relationships and new insights on ancient intergeneric hybridization.

After the merger of the former Taxodiaceae and Cupressaceae s.s., currently the conifer family Cupressaceae (sensu lato) comprises seven subfamilies and 32 genera, most of which are important components of temperate and mountainous forests. With the exception of a recently published genus-level phylogeny of gymnosperms inferred from sequence analysis of 790 orthologs, previous phylogenetic studies of Cupressaceae were based mainly on morphological characters or a few molecular markers, and did not completely resolve the intergeneric relationships. In this study, we reconstructed a robust and well-resolved phylogeny of Cupressaceae represented by all 32 genera, using 1944 genes (Orthogroups) generated from transcriptome sequencing. Reticulate evolution analyses detected a possible ancient hybridization that occurred between ancestors of two subclades of Cupressoideae, including Microbiota-Platycladus-Tetraclinis (MPT) and Juniperus-Cupressus-Hesperocyparis-Callitropsis-Xanthocyparis (JCHCX), although both concatenation and coalescent trees are highly supported. Moreover, divergence time estimation and ancestral area reconstruction indicate that Cupressaceae very likely originated in Asia in the Triassic, and geographic isolation caused by continental separation drove the vicariant evolution of the two subfamilies Cupressoideae and Callitroideae in the northern and southern hemispheres, respectively. Evolutionary analyses of some morphological characters suggest that helically arranged linear-acicular leaves and imbricate bract-scale complexes represent ancestral states, and the shift from linear-acicular leaves to scale-like leaves was associated with the shift from helical to decussate arrangement. Our study sheds new light on phylogeny and evolutionary history of Cupressaceae, and strongly suggests that both dichotomous phylogenetic and reticulate evolution analyses be conducted in phylogenomic studies.

Phylotranscriptomics reveals the evolutionary history of subtropical East Asian white pines: further insights into gymnosperm diversification.

Floristic composition within a geographic area is driven by a wide array of factors from local biotic interactions to biogeographical processes. Subtropical East Asia is a key biodiversity hotspot of the world, and harbors the most families of extant gymnosperms and a large number of endemic genera with ancient origins, but rare phylogenetic studies explored whether it served as a diversification center for gymnosperms. Here, we investigated the evolutionary and biogeographical history of subtropical East Asian white pines using an integrative approach that combines phylotranscriptomic and ecological analyses. Using 2,606 orthologous nuclear genes, we reconstructed a fully resolved and dated phylogeny of these species. Two main clades first diverged in the early Miocene, and by the late Miocene, all species appeared. Two white pines endemic to Taiwan Island experienced independent colonization events and regional extinction, which resulted in the present disjunctive distribution from mainland China. Ecological and biogeographical analyses indicate that the monsoon-driven assembly of evergreen broadleaved forests (EBLFs) might have significantly affected the diversification of subtropical East Asian white pines. Our study highlights the interactions of biotic and abiotic forces in the diversification and speciation of subtropical East Asian white pines. These findings indicate that subtropical East Asia is not only a floristic museum, but also a diversification center for gymnosperms. Our study also demonstrates the importance of phylotranscriptomics on species delimitation and biodiversity conservation, particularly for closely related species.

New insights into the phylogeny and evolution of Podocarpaceae inferred from transcriptomic data.

Phylogenies of an increasing number of taxa have been resolved with the development of phylogenomics. However, the intergeneric relationships of Podocarpaceae, the second largest family of conifers comprising 19 genera and approximately 187 species mainly distributed in the Southern Hemisphere, have not been well disentangled in previous studies, even when genome-scale data sets were used. Here we used 993 nuclear orthologous groups (OGs) and 54 chloroplast OGs (genes), which were generated from 47 transcriptomes of Podocarpaceae and its sister group Araucariaceae, to reconstruct the phylogeny of Podocarpaceae. Our study completely resolved the intergeneric relationships of Podocarpaceae represented by all extant genera and revealed that topological conflicts among phylogenetic trees could be attributed to synonymous substitutions. Moreover, we found that two morphological traits, fleshy seed cones and flattened leaves, might be important for Podocarpaceae to adapt to angiosperm-dominated forests and thus could have promoted its species diversification. In addition, our results indicate that Podocarpaceae originated in Gondwana in the late Triassic and both vicariance and dispersal have contributed to its current biogeographic patterns. Our study provides the first robust transcriptome-based phylogeny of Podocarpaceae, an evolutionary framework important for future studies of this family.

Spatiotemporal Evolution of the Global Species Diversity of Rhododendron.

Evolutionary radiation is a widely recognized mode of species diversification, but its underlying mechanisms have not been unambiguously resolved for species-rich cosmopolitan plant genera. In particular, it remains largely unknown how biological and environmental factors have jointly driven its occurrence in specific regions. Here, we use Rhododendron, the largest genus of woody plants in the Northern Hemisphere, to investigate how geographic and climatic factors, as well as functional traits, worked together to trigger plant evolutionary radiations and shape the global patterns of species richness based on a solid species phylogeny. Using 3,437 orthologous nuclear genes, we reconstructed the first highly supported and dated backbone phylogeny of Rhododendron comprising 200 species that represent all subgenera, sections, and nearly all multispecies subsections, and found that most extant species originated by evolutionary radiations when the genus migrated southward from circumboreal areas to tropical/subtropical mountains, showing rapid increases of both net diversification rate and evolutionary rate of environmental factors in the Miocene. We also found that the geographically uneven diversification of Rhododendron led to a much higher diversity in Asia than in other continents, which was mainly driven by two environmental variables, that is, elevation range and annual precipitation, and were further strengthened by the adaptation of leaf functional traits. Our study provides a good example of integrating phylogenomic and ecological analyses in deciphering the mechanisms of plant evolutionary radiations, and sheds new light on how the intensification of the Asian monsoon has driven evolutionary radiations in large plant genera of the Himalaya-Hengduan Mountains.

NS1: A Key Protein in the "Game" Between Influenza A Virus and Host in Innate Immunity.

Since the influenza pandemic occurred in 1918, people have recognized the perniciousness of this virus. It can cause mild to severe infections in animals and humans worldwide, with extremely high morbidity and mortality. Since the first day of human discovery of it, the "game" between the influenza virus and the host has never stopped. NS1 protein is the key protein of the influenza virus against host innate immunity. The interaction between viruses and organisms is a complex and dynamic process, in which they restrict each other, but retain their own advantages. In this review, we start by introducing the structure and biological characteristics of NS1, and then investigate the factors that affect pathogenicity of influenza which determined by NS1. In order to uncover the importance of NS1, we analyze the interaction of NS1 protein with interferon system in innate immunity and the molecular mechanism of host antagonism to NS1 protein, highlight the unique biological function of NS1 protein in cell cycle.

Both Conifer II and Gnetales are characterized by a high frequency of ancient mitochondrial gene transfer to the nuclear genome.

BACKGROUND: Mitochondrial gene transfer/loss is common in land plants, and therefore the fate of missing mitochondrial genes has attracted more and more attention. The gene content of gymnosperm mitochondria varies greatly, supplying a system for studying the evolutionary fate of missing mitochondrial genes. RESULTS: Here, we studied the tempo and pattern of mitochondrial gene transfer/loss in gymnosperms represented by all 13 families, using high-throughput sequencing of both DNA and cDNA. All 41 mitochondrial protein-coding genes were found in cycads, Ginkgo and Pinaceae, whereas multiple mitochondrial genes were absent in Conifer II and Gnetales. In Conifer II, gene transfer from mitochondria to the nucleus followed by loss of the mitochondrial copy was common, but complete loss of a gene in both mitochondrial and nuclear genomes was rare. In contrast, both gene transfer and loss were commonly found in Gnetales. Notably, in Conifer II and Gnetales, the same five mitochondrial genes were transferred to the nuclear genome, and these gene transfer events occurred, respectively, in ancestors of the two lineages. A two-step transfer mechanism (retroprocessing and subsequent DNA-mediated gene transfer) may be responsible for mitochondrial gene transfer in Conifer II and Gnetales. Moreover, the mitochondrial gene content variation is correlated with gene length, GC content, hydrophobicity, and nucleotide substitution rates in land plants. CONCLUSIONS: This study reveals a complete evolutionary scenario for variations of mitochondrial gene transferring in gymnosperms, and the factors responsible for mitochondrial gene content variation in land plants.

Phylogenomic and ecological analyses reveal the spatiotemporal evolution of global pines.

How coniferous forests evolved in the Northern Hemisphere remains largely unknown. Unlike most groups of organisms that generally follow a latitudinal diversity gradient, most conifer species in the Northern Hemisphere are distributed in mountainous areas at middle latitudes. It is of great interest to know whether the midlatitude region has been an evolutionary cradle or museum for conifers and how evolutionary and ecological factors have driven their spatiotemporal evolution. Here, we investigated the macroevolution of Pinus, the largest conifer genus and characteristic of northern temperate coniferous forests, based on nearly complete species sampling. Using 1,662 genes from transcriptome sequences, we reconstructed a robust species phylogeny and reestimated divergence times of global pines. We found that ∼90% of extant pine species originated in the Miocene in sharp contrast to the ancient origin of Pinus, indicating a Neogene rediversification. Surprisingly, species at middle latitudes are much older than those at other latitudes. This finding, coupled with net diversification rate analysis, indicates that the midlatitude region has provided an evolutionary museum for global pines. Analyses of 31 environmental variables, together with a comparison of evolutionary rates of niche and phenotypic traits with a net diversification rate, found that topography played a primary role in pine diversification, and the aridity index was decisive for the niche rate shift. Moreover, fire has forced diversification and adaptive evolution of Pinus Our study highlights the importance of integrating phylogenomic and ecological approaches to address evolution of biological groups at the global scale.

Phylotranscriptomics reveals the complex evolutionary and biogeographic history of the genus Tsuga with an East Asian-North American disjunct distribution.

The disjunct distribution between East Asia and North America is one of the best established biogeographic patterns. A robust phylogeny is fundamental for understanding the biogeographic histories of taxa with this distribution pattern. Tsuga (hemlock) is a genus of Pinaceae with a typical intercontinental disjunct distribution in East Asia and eastern and western North America, and its phylogeny has not been completely reconstructed in previous studies. In this study, we reconstructed a highly resolved phylogeny of Tsuga using 881 nuclear genes, 60 chloroplast genes and 23 mitochondrial genes and explored its biogeographic and reticulate evolutionary history. The results of phylogenetic analysis, molecular dating and ancestral area reconstruction indicate that Tsuga very likely originated from North America in the late Oligocene and dispersed from America to East Asia via the Bering Land Bridge during the middle Miocene. In particular, we found complex reticulate evolutionary pattern among the East Asian hemlock species. T. sieboldii possibly originated from hybridization with the ancestor of T. chinensis from mainland China and T. forrestii as the paternal donor and the ancestor of T. diversifolia and T. ulleungensis as the maternal donor. T. chinensis (Taiwan) could have originated by hybridization together with T. sieboldii and then evolved independently after dispersal to the Taiwan Island, subsequently experiencing mitochondrial DNA introgression with T. chinensis from mainland China. Moreover, our study found that T. chinensis from western China is more closely related to T. forrestii than to T. chinensis from eastern China. The nonmonophyletic T. chinensis needs taxonomic reconsideration.

Unbiased Subgenome Evolution in Allotetraploid Species of Ephedra and Its Implications for the Evolution of Large Genomes in Gymnosperms.

The evolutionary dynamics of polyploid genomes and consequences of polyploidy have been studied extensively in angiosperms but very rarely in gymnosperms. The gymnospermous genus Ephedra is characterized by a high frequency of polyploidy, and thus provides an ideal system to investigate the evolutionary mode of allopolyploid genomes and test whether subgenome dominance has occurred in gymnosperms. Here, we sequenced transcriptomes of two allotetraploid species of Ephedra and their putative diploid progenitors, identified expressed homeologs, and analyzed alternative splicing and homeolog expression based on PacBio Iso-Seq and Illumina RNA-seq data. We found that the two subgenomes of the allotetraploids had similar numbers of expressed homeologs, similar percentages of homeologs with dominant expression, and approximately equal numbers of isoforms with alternative splicing, showing an unbiased subgenome evolution as in a few polyploid angiosperms, with a divergence of the two subgenomes at ∼8 Ma. In addition, the nuclear DNA content of the allotetraploid species is almost equal to the sum of two putative progenitors, suggesting limited genome restructuring after allotetraploid speciation. The allopolyploid species of Ephedra might have undergone slow diploidization, and the unbiased subgenome evolution implies that the formation of large genomes in gymnosperms could be attributed to even and slow fractionation following polyploidization.

The flattened and needlelike leaves of the pine family (Pinaceae) share a conserved genetic network for adaxial-abaxial polarity but have diverged for photosynthetic adaptation.

BACKGROUND: Leaves have highly diverse morphologies. However, with an evolutionary history of approximately 200 million years, leaves of the pine family are relatively monotonous and often collectively called "needles", although they vary in length, width and cross-section shapes. It would be of great interest to determine whether Pinaceae leaves share similar morpho-physiological features and even consistent developmental and adaptive mechanisms. RESULTS: Based on a detailed morpho-anatomical study of leaves from all 11 Pinaceae genera, we particularly investigated the expression patterns of adaxial-abaxial polarity genes in two types of leaves (needlelike and flattened) and compared their photosynthetic capacities. We found that the two types of leaves share conserved spatial patterning of vasculatures and genetic networks for adaxial-abaxial polarity, although they display different anatomical structures in the mesophyll tissue differentiation and distribution direction. In addition, the species with needlelike leaves exhibited better photosynthetic capacity than the species with flattened leaves. CONCLUSIONS: Our study provides the first evidence for the existence of a conserved genetic module controlling adaxial-abaxial polarity in the development of different Pinaceae leaves.

The complete mitochondrial genome of Taxus cuspidata (Taxaceae): eight protein-coding genes have transferred to the nuclear genome.

BACKGROUND: Gymnosperms represent five of the six lineages of seed plants. However, most sequenced plant mitochondrial genomes (mitogenomes) have been generated for angiosperms, whereas mitogenomic sequences have been generated for only six gymnosperms. In particular, complete mitogenomes are available for all major seed plant lineages except Conifer II (non-Pinaceae conifers or Cupressophyta), an important lineage including six families, which impedes a comprehensive understanding of the mitogenomic diversity and evolution in gymnosperms. RESULTS: Here, we report the complete mitogenome of Taxus cuspidata in Conifer II. In comparison with previously released gymnosperm mitogenomes, we found that the mitogenomes of Taxus and Welwitschia have lost many genes individually, whereas all genes were identified in the mitogenomes of Cycas, Ginkgo and Pinaceae. Multiple tRNA genes and introns also have been lost in some lineages of gymnosperms, similar to the pattern observed in angiosperms. In general, gene clusters could be less conserved in gymnosperms than in angiosperms. Moreover, fewer RNA editing sites were identified in the Taxus and Welwitschia mitogenomes than in other mitogenomes, which could be correlated with fewer introns and frequent gene losses in these two species. CONCLUSIONS: We have sequenced the Taxus cuspidata mitogenome, and compared it with mitogenomes from the other four gymnosperm lineages. The results revealed the diversity in size, structure, gene and intron contents, foreign sequences, and mutation rates of gymnosperm mitogenomes, which are different from angiosperm mitogenomes.

Phylogenomics disentangles the evolutionary history of spruces (Picea) in the Qinghai-Tibetan Plateau: Implications for the design of population genetic studies and species delimitation of conifers.

A laborious and difficult task in current tree of life reconstruction is to resolve evolutionary relationships of closely related congeneric species that originated from recent radiations. This is particularly difficult for forest species with long generation times and large effective population sizes such as conifers. The Qinghai-Tibetan Plateau (QTP) and adjacent areas are considered a species diversity center of Picea, harboring 11 species (including 5 varieties) of this genus, but evolutionary relationships of these species are far from being resolved due to recent radiations, morphological convergence, and frequent interspecific gene flow. In this study, we use these spruce species to test whether phylotranscriptomic analysis, combined with population genetic analysis, can disentangle their evolutionary relationships, and to explore whether reticulate evolution has occurred among them. Phylogenomic analyses indicate that all spruce species in the QTP and neighboring areas, except P. asperata and P. crassifolia, cluster together, and in particular, nearly all taxa (including varieties) reflect reciprocally monophyletic lineages, although the two species P. likiangensis and P. brachytyla are not monophyletic. We found that, compared to herbaceous plants, many more genes (a minimum of 600 OGs for Picea) are required to resolve interspecific relationships of conifers. Contrary to previous studies, our data do not support a hybrid origin of P. purpurea, but suggests a hybrid origin for P. brachytyla var. brachytyla and P. likiangensis var. rubescens. We emphasize that the species or species complex used for population genetic and phylogeographical studies should be monophyletic.

Phylotranscriptomics resolves interspecific relationships and indicates multiple historical out-of-North America dispersals through the Bering Land Bridge for the genus Picea (Pinaceae).

A robust phylogeny is prerequisite to understand the evolution and biogeography of organisms. However, ancient and recent evolutionary radiations occurred in many plant lineages, which pose great challenges for phylogenetic analysis, especially for conifers characterized by large effective population sizes and long generation times. Picea is an important component of the dark coniferous forests in the Northern Hemisphere. Previous studies improved our understanding of its evolutionary history, but its interspecific relationships and biogeographic history remain largely unresolved. In the present study, we reconstructed a well-resolved phylogeny of Picea by comparative transcriptomic analysis based on a complete species sampling. The phylogenetic analysis, together with molecular dating and ancestral area reconstruction, further supports the North American origin hypothesis for Picea, and indicates that this genus experienced multiple out-of-North America dispersals by the Bering Land Bridge. We also found that spruces in the Japanese Archipelago have multiple origins, and P. morrisonicola from the Taiwan Island has a close relationship with species from the Qinghai-Tibetan Plateau and adjacent regions. Our study provides the first complete phylogeny of Picea at the genomic level, which is important for future studies of this genus.

Cryptic speciation in the Chinese white pine (Pinus armandii): Implications for the high species diversity of conifers in the Hengduan Mountains, a global biodiversity hotspot.

Conifers are the largest and ecologically and economically most important component group of the gymnosperms. Despite their slow rate of molecular evolution, rapid and recent diversification was unexpectedly prevalent in this ancient group in the Hengduan Mountains, a world's biodiversity hotspot and gymnosperm diversity center in Southwest China. In this study, we investigated the underlying mechanisms and disentangled the interactions of geography and ecology in speciation and evolution in Pinus armandii, an important forest tree species endemic to China, by integrating analyses of population transcriptomics, population genetics and ecological niche modeling. Many lines of evidence suggest that cryptic speciation has occurred in P. armandii. During the process, geologically induced formation of Mount Gongga and other massive peaks might trigger the initial vicariance isolation of the northern and southern subdivisions, and ecologically based selection then reinforced their differentiation and local adaptation. Our ecological niche analysis and earlier reciprocal transplant experiments in P. armandii provided convincing evidences for the critical role of ecology in the process of speciation. These findings suggest that both geography and ecology contributed significantly to the abundance of very recent and rapid species divergences, which promoted the rising of the extremely high conifer diversity in the Hengduan Mountains.

Phylogeny and evolutionary history of Pinaceae updated by transcriptomic analysis.

Pinaceae comprises 11 genera, and represents the largest family of conifers with an extensive wild distribution in the Northern Hemisphere. Intergeneric relationships of Pinaceae have been investigated using many morphological characters and molecular markers, but phylogenetic positions of four genera, including Cathaya, Cedrus, Nothotsuga and Pseudolarix, remain controversial or have not been completely resolved. To completely resolve the intergeneric relationships of Pinaceae, we conducted a comparative transcriptomic study of 14 species representing all Pinaceae genera. Multiple data sets, containing up to 6,369,681 sites across 4676 loci, were analyzed using concatenation and coalescent methods. Our study generated a robust topology, which divides Pinaceae into two clades, one (pinoid) including Cathaya, Larix, Picea, Pinus, and Pseudotsuga, and the other (abietoid) including Abies, Cedrus, Keteleeria, Nothotsuga, Pseudolarix, and Tsuga. Cathaya and Pinus form a clade sister to Picea; Cedrus is sister to the remaining abietoid genera, and the two genera Nothotsuga and Tsuga form a clade sister to Pseudolarix. The discordant positions of Cathaya, Cedrus and Pseudolarix in different gene trees could be explained by ancient radiation and/or molecular homoplastic evolution. The hybrid origin hypothesis of Nothotsuga is not supported. Based on molecular dating, extant Pinaceae genera diverged since about 206 Mya, earlier than the break-up of Pangea, and the divergence among the pinoid genera occurred earlier than the split among the abietoid genera. Moreover, our study indicates that two radiation events occurred in the evolution of Pinaceae genera, and some important morphological characters evolved multiple times based on ancestral state reconstruction.

Phylogenomics resolves the deep phylogeny of seed plants and indicates partial convergent or homoplastic evolution between Gnetales and angiosperms.

After decades of molecular phylogenetic studies, the deep phylogeny of gymnosperms has not been resolved, and the phylogenetic placement of Gnetales remains one of the most controversial issues in seed plant evolution. To resolve the deep phylogeny of seed plants and to address the sources of phylogenetic conflict, we conducted a phylotranscriptomic study with a sampling of all 13 families of gymnosperms and main lineages of angiosperms. Multiple datasets containing up to 1 296 042 sites across 1308 loci were analysed, using concatenation and coalescence approaches. Our study generated a consistent and well-resolved phylogeny of seed plants, which places Gnetales as sister to Pinaceae and thus supports the Gnepine hypothesis. Cycads plus Ginkgo is sister to the remaining gymnosperms. We also found that Gnetales and angiosperms have similar molecular evolutionary rates, which are much higher than those of other gymnosperms. This implies that Gnetales and angiosperms might have experienced similar selective pressures in evolutionary histories. Convergent molecular evolution or homoplasy is partially responsible for the phylogenetic conflicts in seed plants. Our study provides a robustly reconstructed backbone phylogeny that is important for future molecular and morphological studies of seed plants, in particular gymnosperms, in the light of evolution.